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; ABSTRACT 

"iJ . We present the 15 /urn extragalactic source counts from the Final Analysis Catalogue 

' of the European Large Area ISO Survey southern hemisphere field SI, extracted using 

\ the Lari method. The large number of extragalactic sources (~ 350) detected over this 

— . i area between about 0.5 and 100 mJy guarantee a high statistical significance of the 

' source counts in the previously poorly covered flux density range between IRAS and 

the Deep ISOCAM Surveys. The bright counts in SI (£15 ^ > 2 mJy) are significantly 
lower than other published ISOCAM counts in the same flux range and are consistent 
S-^ I with a flat, Euclidean slope, suggesting the dominance of a non-evolving population. 

c/3 . In contrast, at fainter fluxes (S15 < 2 mJy) our counts show a strong departure 

from no evolution models, with a very steep super-Euclidean slope down to our flux 
limit (^0.5 mJy). Strong luminosity and density evolution of the order of, respectively, 
L oc (1 + z) 3 and p cx (1 + z) 3 5 is needed at least for the population of star-forming 
galaxies in order to fit the counts and the redshift distributions observed at different 



fluxes. A luminosity break around 10 10 ' 8 Lq must be introduced in the local luminosity 
function of starburst galaxies in order to reproduce our sharp increase of the counts 
below 2 mJy and the redshift distributions observed for 15 /^m sources at different 
flux levels. The contribution of the strongly evolving starburst population (down to 
50 pJy) to the 15 /xm cosmic background is estimated to be ~2.2 nW m -2 sr _1 , which 
corresponds to ~67% of the total mid-infrared background estimate. 

Key words: galaxies: evolution - galaxies: starburst - cosmology: observations- 
infrared: galaxies. 



1 INTRODUCTION 

Deep galaxy counts are a key instrument for the study of 
galaxy evolution and can provide strong constraints to the- 
oretical models. In fact, the departure of source counts from 
Euclidean predictions depends on the intrinsic evolution of 
galaxies and their redshift distribution. In the past few years 
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several deep observations in different wave-bands have pro- 
vided a significant advance in our knowledge of galaxy for- 
mation and evolution. Optical surveys found a strong evolu- 
tion of the population of blue galaxies as a function of red- 
shift to z ~ 1 (Lilly et al. 1996; Metcalfe et al. 1995), while 
Mid- and Far-Infrared deep surveys (Elbaz et al. 1999a,b, 
2000; Dole et al. 1999), together with the detection of a 
substantial diffuse cosmic Infrared Background in the 300 
fim - 1 mm range (Puget et al. 1996; Hauser et al. 1998; 
Fixsen et al. 1998; Lagache et al. 1999) implied a strong evo- 
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lution also for galaxies emitting in the infrared. In fact, the 
Mid/Far-IR extragalactic background is at least as large as 
the UV/optical/NIR background, thus implying a stronger 
contribution of obscured star formation at redshifts larger 
than those observed by IRAS. 

IRAS has sampled the local Universe (z < 0.2) in the 
Mid/Far-IR band, discovering Luminous Infrared Galax- 
ies (LIGs: L > 10 11 1/© ) which radiate most of their light 
in the infrared band. Although LIGs are the most lumi- 
nous starburst galaxies ever detected, they are relatively 
rare in the local Universe, thus making up only a small 
fraction of the total energy output from galaxies (Soifer & 
Negeubauer 1991). However, different analyses of the IRAS 
extragalactic source counts by Hacking, Houck and Condon 
(1987), Lonsdale & Hacking (1989), Saunders et al. (1990) 
and Kim and Sanders (1998) have shown some evidence for 
strong evolution at low flux density levels for ULIGs (Ultra- 
LIGs: LboI — Lir > 1Q 12 Lq). Due to the small redshift 
range sampled by IRAS, these results can only be indica- 
tive, though the suggestion that ULIGs might have played 
a stronger role in the past is supported by the detection of 
the strong infrared background. 

With a thousand times better sensitivity and sixty times 
better resolution than IRAS, ISOCAM instrument (Ce- 
sarsky et al. 1996) on board of the Infrared Space Obser- 
vatory (ISO; Kessler et al. 1996) has provided Deep and 
Ultra-Deep mid-infrared extragalactic surveys (mainly with 
the LW3 filter: 12 - 18 pm), unveiling most of the star- 
formation in the Universe to z — 1. The source counts 
derived from these Deep/Ultra-Deep surveys (covering the 
flux density range 0.05 - 4 mjy) strongly diverge from no- 
evolution models at fluxes fainter than about 1 mjy, with an 
increasing difference that reaches a factor of 10 around the 
faintest limits (0.05 - 0.1 mjy; Elbaz et al. 1999a,b). The 
faint mid- infrared sources detected in the Deep/Ultra-Deep 
ISOCAM surveys have been identified mainly with galax- 
ies at z ~ 0.7 and show LIG-like luminosities (Aussel et al. 
1999b; Elbaz et al. 1999a,b; Flores et al. 1999). 

Although the Deep/Ultra-Deep ISOCAM Surveys have 
produced crucial results on galaxy evolution in the infrared, 
having identified most of the galaxies producing the mid- 
infrared background, there is a large gap in the flux density 
sampled by these surveys and by the IRAS Surveys. In par- 
ticular, the flux density range 4 - 200 mjy is essentially un- 
covered in the Mid/Far-infrared and only a few sources have 
been detected by the Deep ISOCAM Surveys in the impor- 
tant flux range 1-4 mjy, where most of the existing evolu- 
tionary models (i.e. Franceschini et al. 2001; Xu 2000; Chary 
& Elbaz 2001; Rowan-Robinson 2001) predict a substantial 
change in the relative contribution of a local non-evolving 
and a more distant evolving population. This scarcity of 
data in this flux interval reflects also in a poor knowledge 
of the local luminosity function for starburst galaxies and a 
corresponding uncertainty in the evolutionary properties of 
the different classes of infrared extragalactic objects. 

The European Large Area ISO Survey (ELAIS; Oliver 
et al. 2000) is the largest survey conducted with the Infrared 
Space Observatory (ISO) and provides a link between the 
IRAS survey and the Deep/Ultra-Deep ISOCAM surveys. 
ELAIS is a collaboration between 20 European institutes 
which involves a deep, wide-angle survey at high galactic 
latitudes, at wavelengths of 6.7 (im (LW2), 15 fim (LW3), 



90 jum (C100) and 175 ^m (C200) with ISO. In particular, 
the 15 fim survey covers a total area of ~ 12 deg 2 , divided 
into 4 main fields and several smaller areas. One of the main 
fields, SI, is located in the southern hemisphere. The whole 
SI area has been surveyed in the radio (at 1.4 GHz, Grup- 
pioni et al. 1999) , in several optical bands (La Franca et al. 
2002, in preparation) and in the hard X-ray with BeppoSAX 
(Alexander et al. 2001). Moreover, spectroscopic informa- 
tion and redshifts are available for a large number of sources 
(Gruppioni et al. 2001; La Franca, Gruppioni, Matute et al. 
2002, in preparation). 

We have reduced the 15^tm data in SI using a new ISO- 
CAM data reduction technique (LARI technique) especially 
developed for the detection of faint sources, obtaining a cat- 
alogue (complete at the 5a level) of 462 sources in the flux 
density range 0.45 - 150 mjy. Details about the data reduc- 
tion technique and the source catalogue have been presented 
in Lari et al. 2001 (hereafter Paper I). Here we present the 
source counts at 15 /im derived from that catalogue. The 
paper is structured as follows. In section ^ we give a brief 
description of the 15 /im ELAIS survey in SI. In section ^ 
we present the completeness and reliability of our sample at 
different flux levels. In section ^ we present the sample used 
to derive the source counts, which are shown in section ^| 
Finally, in sections ^ and |^ we discuss our results and their 
implications and present our conclusions. 

Throughout this paper we will assume Ho = 50 km s _1 
Mpc" 1 , Q, m = 0.3 and fi A = 0.7. 



2 DESCRIPTION OF THE ELAIS SI SURVEY 

The ELAIS survey at 15 /im, performed in raster mode with 
the ISOCAM instrument on board of ISO, covers a total area 
of ~ 12 deg 2 divided into 4 main fields and several smaller 
areas. The main field located in the southern hemisphere 
(SI) is centered at a(2000) = 00 h 34 m 44.4 s , (5(2000) = 
-43° 28' 12" and covers an area of 2° x 2° . The 15 /im survey 
performed in SI with the ISOCAM instrument consists of 9 
different rasters. Each raster covers an area of ~ 43.5' x 42'; 
eight of them have been observed once, while one, Sl_5, was 
observed three times. 

The 15 £tm data have been reduced and analyzed using 
the LARI technique, described in detail in Paper I. With 
this data reduction method, we have obtained a sample of 
462 sources with signal-to-noise ratio > 5 in the flux range 
0.45 - 150 mjy. The fainter sources have been detected in 
the central raster of SI (Sl_5), whose image has been ob- 
tained by combining three single observations centered on 
the same position. The source catalogue in SI and the rela- 
tive parameter errors have been presented and discussed in 
Paper I, together with the detection rates at different flux 
densities derived with simulations. 

The detection rates given in Paper I cannot be directly 
translated to completeness of the real catalogue, because our 
simulations were performed at discrete flux values rather 
than following a continuous flux distribution. However, the 
results presented in Paper I can be used to obtain the com- 
pleteness of the catalogue and the source counts corrections, 
as discussed in the next section. 
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3 COMPLETENESS AND RELIABILITY 

3.1 Brief Summary of the Results and Definitions 
from Paper I 

Before describing in detail the method used to derive the 
completeness of our source counts, it is useful to summa- 
rize the more relevant results and definitions of Paper I. 
The simulations performed in the SI field provided not only 
the completeness of our detections at different flux levels, 
but also the internal calibration of the source photometry 
and the distribution of the ratio between the measured and 
the theoretical peak fluxes (crucial for the computations de- 
scribed in next section). 

Here we give some relevant definitions and relations: 

• f 3 : is the peak flux measured on maps for both real and 
simulated sources. Its value depends both on data reduc- 
tion method and on ELAIS observing strategy + ISOCAM 
instrumental effects; 

• fo' is the 'theoretical' peak flux measured on simulated 
maps containing neither glitches or noise. Its value depends 
only on ELAIS observing strategy + ISOCAM instrumental 
effects; 

• ? = fs/h; 

• q-med'- is the peak of the f 3 /fo distribution (also called 
systematic flux bias) and is 0.78 ± 0.03 in SI and 0.82 ± 
0.03 in Sl_5. These values are used to correct the measured 
flux densities; 

• flux density determination: 

S = \Jjr X s ( ^j /q m ed (1) 

where for simulations so is the injected total flux, while for 
real data so is derived through successive iterations starting 
from a rough estimated value: 

so = fs/{f a /s } sim (2) 

where {fa/so) sim = 0.132 is the average value resulting from 
simulations. 

We can consider s as the measured flux density and so as 
the 'true' flux density of a source; 

• the photometric accuracy of our data reduction of the 
SI area has been tested using the stars of the field and fol- 
lowing the relation calibrated on IRAS data by Aussel & 
Alexander (2001) to predict the fluxes of stars. By a com- 
parison between the predicted fluxes and the ones derived 
from our analysis, we have obtained a very good agreement 
and a relative flux scale of 1.096 ± 0.044 (i.e. our fluxes 
have to be multiplied by 1.096 to be put on the same scale 
as IRAS fluxes). 

3.2 The g Function 

As described in Paper I and mentioned in the previous sec- 
tion, with simulations in SI we have derived the distribution 
of measured (f a ) to theoretical (/o) peak flux ratio. This dis- 
tribution is crucial in deriving the completeness of the cata- 
logue and the internal flux calibration. To this purpose, we 
have considered the / s //o distribution as a model function, 
hereafter called g function. The g distribution function, ob- 
tained for simulated sources, is a combination of an intrinsic 
g function (referred as go) plus a term due to noise. As a 



rough estimate of go in Paper I we considered the measured 
g distribution obtained for 3 mjy simulated sources (<7o-3)- 
Because 3 mjy is a relatively high flux density (it is the 
highest flux injected in our simulations), the gos distribu- 
tion is almost unbiased by detection incompleteness, and we 
can consider it to be a good approximation of the intrinsic 
one. However, using the same simulation, we have obtained a 
more refined estimate of the intrinsic go function correcting 
the observed fs/fo distribution for the low level of detection 
incompleteness still present at this flux. Our procedure was 
the following: 

each detected source i has a c/i = f 3 (i) / fo (i) value. On a dif- 
ferent position j, assuming that the fs(i)/ fo(i) value is not 
depending on position, the peak flux would be measured 
with a value f 3 ,i(j) = qi x fo(j) and would be detected only 
if f s .i{j) were greater than 5 times the local noise (c(j)). If 
Ni is the number of possible detections in all the different 
j positions of the simulations, the weight to be assigned to 
the i th source is l/Ni. Finally, the estimate of the intrinsic 
go function in a given interval of q is obtained by summing 
the weights of all the detected sources with a q value in the 
same interval. By construction, the integral of go over the 
entire range of q is one. 

Once obtained the go distribution, the general g distribution 
can be derived by convolving go with the distribution of 
noise, here assumed to be Gaussian. This is possible in the 
approximation that the measured peak flux f s is the sum of 
a "true" value, f s , and a stochastic term, \, due to noise, so 
that the q (= fs/fo) measured value is: 

fa fs | X /o\ 

q = -r = T + T ( 3 ) 
Jo Jo fo 

with the x term following a Gaussian distribution with av- 
erage equal to zero and dispersion equal to the rms noise 
value, o. Following this formalism, then the g distribution, 
expressed as a function of fs/fo and fo/o, is obtained by 
just convolving go with the noise distribution, integrating 
over the possible values of the variable \- 

^) = Jl~ *X (4) 

fo cr V 2tt cr J f fo 

In figure |l| we show the intrinsic go distribution computed 
as above (solid histogram) and the predicted distributions 
of the ratio f s /fo in presence of a a of 26 /ijy/pixel (typical 
of our data) for three different values of fo, corresponding 
to different mean values of fo for different total input fluxes. 
As we can notice, the presence of noise broadens the flux dis- 
tributions and this effect becomes stronger towards fainter 
fluxes (as shown also in figure 7 of Paper I) . 
If we assume that the g function reflects all the multiplica- 
tive and additive error components due to the data reduc- 
tion, we can use this function to predict the distribution of 
the detection rate for both simulated and real sources, as 
described in the next subsection. 



3.3 Completeness 

First we have computed the incompleteness introduced by 
our data reduction method, represented by the loss of 
sources not accounted for by our model. In fact, sources can 
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Figure 1. Predicted distributions of the ratio between the mea- 
sured and the theoretical peak flux, f s /fo (also g function), in 
presence of a rms noise of 26 /zjy/pixel. These distributions, ob- 
tained from simulations, are shown for three different average 
values of fo corresponding to different input fluxes. The solid his- 
togram shows the intrinsic go distribution (see text) 



be missed by our method if interpreted as background tran- 
sients. The incompleteness of our method is obtained from 
simulations by computing the ratio between the number of 
detections and the number of expected sources in different 
peak flux intervals ([fo — dfo, fo + dfo]). The number of ex- 
pected sources is derived by summing together the predicted 
detection contributions (h(fo) = L ,, g(q, fo)dq) of all the 

J OCT/ Jo 

sources i with peak flux fo% belonging to the same interval: 

N exp (f ) = /l(/o) (5) 

/o.e[/o±d/ ] 

In figure |i| the resulting function describing the incom- 
pleteness of our method is plotted as a function of fo, to- 
gether with its lower and upper envelopes, for SI (top panel) 
and Sl_5 (bottom panel). 

To obtain the global correction to be applied to our 
source counts, we need to consider also the areal coverage of 
our survey (i.e. the fraction of the survey area where a source 
of peak flux f s can be detected: f s > 5<r) and the fact that 
real sources do not follow a discrete flux distribution as our 
simulated sources. To account for the latter effect, we have 
assumed a certain shape to describe the real source counts 
observed in the sky. According to the published counts at 
15/^m (e.g. Elbaz et al. 1999a), we have assumed two power 
laws between 0.4 and 150 mjy: 



S~ ai if S > 2 mjy 
s - a2 if s < 2 mJy 



(6) 



with ai=2.3 and 6*2=3.0 as first estimate values. 
By weighting the above 'theoretical' counts per unit of area 
by the g function convolved with the areal coverage func- 
tion, and with the function describing the completeness 
of our method, we have computed the counts predicted 
in our Survey. Being dN(s) = dN(s)/ds ■ ds the number 
of sources detected in the flux density bin ds, dAf(so) = 
dJ\f(so)/dso ■ dso the theoretical number of sources detected 
in the 'true' flux density bin dso and dso /ds — q m ed /q (from 
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Figure 2. Completeness function of Lari method as function of 
the theoretical peak flux for SI (top) and Sl_5 (bottom). The 
dashed lines give an estimate of the lower and upper envelopes of 
the completeness function. 



equation |l|), we have 
dN 



ds 



(«) = 



dq 



dAf . 
ds 


dso~ 
ds . 


dAf. . 
dso 


qmed 

q 



where gA is the convolution of the g function with the areal 
coverage function and C is the completeness function shown 
in figure ^| The averaging is performed over all the predicted 
peak flux values fo and with a random sampling of source 
positions in the survey. These expected source counts thus 
take into account the effects produced by the specific ob- 
servational parameters of the ELAIS 15 (jm survey and by 
our data reduction method. By putting aVV(so) /dso in the 
form given by equation ^ into equation ^, and changing so 
to s ■ q me d/q (from equation |l|) , we obtain the source counts 
predicted for our survey, dN(s)/ds. The ratio between these 
counts and the theoretical source counts, dAf(so) /dso, gives 
the global correction (including incompleteness and areal 
coverage) to be applied to our measured source counts. The 
inverse of the global correction ( effective area) is reported in 
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Figure 3. Effective area versus total flux density for the 15 /im 
catalogues in SI main area (top) and Sl_5 (bottom). 

figure^] for both SI (top panel) and Sl_5 (bottom panel) sur- 
veys. Since our source counts corrected for incompleteness 
were significantly different from the original model counts 
(especially the power law at faint flux densities), we have 
iterated this procedure by adjusting the power laws param- 
eters until convergence is achieved. The final estimate of the 
model counts dj\f(s)/ds is obtained for «i=2.3 and 02=3.6. 



4 THE DATA SAMPLE USED FOR THE 
SOURCE COUNTS 

The catalogue used to derive the source counts is not ex- 
actly the catalogue published in Paper I, but contains some 
differences, as described in this section. 

First, we have conservatively decided to exclude from 
the source list used for counts 35 sources detected in SI, 
all with flux density < 1.5 mjy, resulting "dubious" at a 
visual inspection of their pixel history. These sources, de- 
tected above the 5<r threshold on the maps obtained through 
a combination of several images, are too faint to be dis- 
tinguished from noise on the single pixel histories without 
uncertainty. Moreover, an additional factor which strength- 
ened our doubts about the reliability of these sources is their 
very low optical identification rate. In fact, while for the en- 



Figure 4. Correction factors as function of signal to noise to be 
applied to the SI (solid) and Sl-5 (dashed) source flux density to 
account for flux loss caused by positional errors in the theoretical 
peak flux fo- 



tire catalogue (minus the 35 "dubious" sources) the optical 
identification rate (within a circle of 4 arcsec radius) on the 
DSS2 images for sources fainter than 1.5 mjy is 60%, for the 
35 "dubious" sources it is only 14%, with a chance detection 
rate of about 10%. 

Second, before computing the source counts we have 
applied three small further corrections to the flux densities 
presented in the catalogue of Paper I: 

1. We have applied the calibration factor of 1.096 derived 
from the comparison with stars in order to put our fluxes on 
the same scale as IRAS fluxes (see section i.l and section 7 
in Paper I). 

2. We have corrected for the average underestimate of 



the true flux introduced by positional errors of the auto- 
simulated peak flux fo (see section 6 in Paper I); Since fo is 
computed on the measured positions and not on the 'true' 
ones, the measured flux is on average underestimated. This 
underestimate is flux dependent, due to the flux dependence 
of positional uncertainties. In figure ^ the corrections to be 
applied to the source flux densities due to this effect are 
reported as a function of signal to noise for both 5*1 (solid 
line) and 5*1_5 (dashed line). 

3. In Sl_5 only we have corrected for the additional loss of 
flux due to the combination of the three rasters. To estimate 
this correction factor, for each source found in the combined 
51_5 map we have measured the flux in the three separate 
rasters and compared their average with the flux measured 
on the combined map. The mean ratio between the combined 
and the averaged single fluxes, considered as the correction 
factor, is R = 0.96 ± 0.01. 

In the following statistical analysis we have chosen not 
to eliminate the 20 repeated sources (those belonging to 
the overlapping regions of two different rasters), but to con- 
sider them as different sources. The reason of this choice is 
that the detectability and completeness analysis, and conse- 
quently the correction factors to be applied to our data, have 
been performed on the single raster areas. This can explain 
why the effective area of SI (see figure ^) corresponding to 
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Figure 5. Magnitude distribution of our stellar identifications Figure 6. ISO-optical separation distribution for the stellar iden- 

(line-filled histogram) to Bj = 16, compared to the magnitude tifications in SI. 
distribution of all the GSC-II stars in the SI area (empty his- 
togram). 



bright fluxes (= 8 times the area of one raster) is larger 
than the area of sky effectively covered by the SI rasters. 
The sum of the effective areas shown in figure ^| is about 4.6 
square degrees, while the SI + Sl_5 survey covers an area 
of 4 square degrees (about 15% of the area is covered by at 
least two rasters). 

Finally, to compute the extragalactic source counts, we 
have excluded from our lists all the sources with a stellar 
counterpart brighter than J = 16 in the Guide Star Cata- 
logue 11^ and with an evident stellar appearance (i.e. point- 
like with spikes) on the DSS2^ images. We have chosen not 
to exclude any stellar identification fainter than Bj — 16 
found in the GSC-II (which is complete to Bj = 19.5), be- 
cause the reddest faint stars in our sample have estimated 
Bj magnitude of the order of 15 (Aussel and Alexander 
2001), and at fainter magnitudes the elimination from the 
sample of stellar-like objects might cause the elimination 
of AGN instead of stars. In fact, at Bj > 16 the identifica- 
tions with point-like objects are expected to be dominated 
by AGN. 

In the end, we have identified with stars, and subtracted 
from our list, 82 sources in SI and 20 in Sl_5 (in total 102 

t The Guide Star Catalogue-II is a joint project of the Space 
Telescope Science Institute and the Osservatorio Astronomico 
di Torino. Space Telescope Science Institute is operated by the 
Association of Universities for Research in Astronomy, for the 
National Aeronautics and Space Administration under contract 
NAS5-26555. The participation of the Osservatorio Astronomico 
di Torino is supported by the Italian Council for Research in 
Astronomy. Additional support is provided by European South- 
ern Observatory, Space Telescope European Coordinating Facil- 
ity, the International GEMINI project and the European Space 
Agency Astrophysics Division. 

t The "Second epoch Survey" (DSS2) of the Southern Sky was 
made by the AAO with the UK Schmidt Telescope. Plates from 
this survey have been digitized and compressed by the STScI. The 
digitized images are copyright(c) 1993-1995 by the AAO Board 
and are distributed herein by agreement. All rights reserved. 
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Figure 7. Ratio between stars and the total number of 15 fim 
sources as a function of flux density for our SI sample. 



stars, 87 of which are different), thus leaving a total of 325 
extragalactic sources in S1+S1_5 (320 different). The Bj 
magnitude (from GSC-II) distribution of these stars is re- 
ported in figure ^ (filled histogram), where also the magni- 
tude distribution of all the GSC-II stars in the SI area is 
reported (to Bj = 16). The maximum separation we find 
between ISO sources and star positions for our stellar iden- 
tifications is < 4 arcsec, as shown in figure ^| Given these 
positional differences and the surface density of stars, we 
estimate that less than one ISO /star association could be 
spurious (down to the considered magnitude limit Bj — 16) . 
In figure M the fraction of stars in SI is plotted as a func- 
tion of flux density. At S15 ^ m > 50 mJy all our sources are 
identified with stars and still at Si 5 Mm — 5 the fraction of 
stars is > 50%. For this reason, an accurate star subtrac- 
tion is very important before computing the extragalactic 
ISOCAM counts, even at relatively faint fluxes (1-2 mJy), 
where the fraction of stars is still of the order of 20%. 

In the end, the extragalactic sample used to compute 
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the source counts is composed by 325 source, 320 of which 
are different. 



5 SOURCE COUNTS 

Given the extent of the 15 /xm survey in SI and the signifi- 
cant depth reached in its central area, Sl_5, our source sam- 
ple is optimally suited to study the ISO CAM source counts 
with a large statistics and over a broad flux range (0.5 - 100 
mjy). Therefore, the combined sample of our Sl(main area) 
+ Sl_5 non-stellar sources with Si5 M m > 5<r has been used to 
construct the mid-infrared extragalactic source counts dis- 
tribution. 

We used the effective areas derived in section | (see fig- 
ure^), to obtain the extragalactic mid-infrared source counts 
from our 15 fim samples. In Table [l] the 15 /im source counts 
in SI, Sl_5 and S1+S1_5 areas are presented. We have first 
computed the counts for the two samples in SI and Sl_5 
separately, as reported in the first eight columns of Table 
[jj, giving respectively the adopted flux density intervals, the 
average flux density in each interval (computed as the ge- 
ometric mean of the two flux density limits), the observed 
number of sources, the effective area (figure ^top), the dif- 
ferential source counts and their associated errors in SI and 
the observed number of sources, the effective area (figure ^ 
bottom) and the differential source counts with their errors 
in Sl_5. The differential source counts have been obtained 
by weighting each single source for its effective area rather 
than weighting the total number of sources in each flux den- 
sity bin for the effective area corresponding to the reference 
flux density of that bin. The errors associated to the counts 
in each bin have been computed as -^/Ei(l / A\ s f{Si)) , where 
the sum is for all the sources with flux density Si belonging 
to the bin and A e ff(Si) is the effective area corresponding 
to that source flux. These errors take into account only the 
Poissonian term of the uncertainties associated to the source 
counts, for consistency with the other literature works. Es- 
pecially at faint flux density, where the effective area (and 
consequently the correction factor) is a steep function of flux, 
the errors quoted in Table [j] should be considered as lower 
limits of the 'true' errors (including also the uncertainty in 
the effective area computation, shown as dashed curves in 
figure |). 

The counts computed for SI and Sl_5 are consistent within 
the errors in most flux density intervals, the only exception 
being the 0.8 - 1.2 mjy flux bin, where the counts in SI 
are larger than the ones in Sl_5 at a formal level of ~ 2.3<r. 
However, since the low flux density errors are somewhat un- 
derestimated, we can consider the source counts in SI and 
Sl_5 consistent in all the common flux density intervals, over 
the whole range 0.5 - 100 mjy. 

For this reason, we have also computed the source counts 
in the "combined" sample (S1+S1_5), by considering all the 
sources as belonging to a unique sample: for each source 
the combined effective area is the sum of the effective ar- 
eas in SI and Sl_5 (whose values are reported in column 
4 and 7 respectively). In the last four columns the total 
number of sources in each flux density bin, the differential 
source counts, the differential counts normalized to the Eu- 
clidean distribution (by multiplying by (s) 2 ' 5 ) with their 
errors and the integral source counts (with errors) for the 



combined sample are reported. In the first flux density bin 
the counts for the combined sample coincide with the Sl-5 
counts, since, due to their negligible effective area, we have 
not considered the SI data. 

Note that the flux bins are partially overlapping, there- 
fore they are not statistically independent (they are alter- 
nately independent). The choice of partially overlapping flux 
density bins for our source counts representation is based on 
the need of a tight sampling of the region where the counts 
start diverging from no evolution models, in order to bet- 
ter determine the break point and the counts shape (with a 
larger statistics). As mentioned in the previous section, in 
computing the counts we have considered as two different 
sources those appearing in two different rasters (belonging 
to the border part of a raster, overlapping with an adjacent 
raster), by suitably weighting them for their detectability 
area in each raster. In fact, in deriving the areal coverage 
function we have considered the total area of each single 
raster, including the overlapping regions. 

The 15 /im differential source counts of the combined 
ELAIS SI and Sl_5 data, normalized to those expected in a 
Euclidean geometry by dividing by S~ 2 ' 5 , are shown in fig- 
ure^ (filled stars). For comparison, source counts from other 
ISOCAM surveys (A2390 from Altieri et al. 1999; ISO HDF- 
N from Aussel et al. 1999a; ISO HDF-S, Marano Firback, 
Marano Ultra-Deep, Marano Deep, Lockman Deep, Lock- 
man Shallow from Elbaz et al. 1999b; data kindly provided 
by D. Fadda, private communication) are also plotted. Our 
counts are lower than the Lockman Deep and Shallow ones. 
However, they appear consistent with the counts obtained 
in the Marano Deep Survey, at least in the common flux 
density range (0.5 - 2 mjy). 

In figure ^ our integral extragalactic source counts are 
reported. The grey-shaded area shows the integral extra- 
galactic counts with 68% confidence contours obtained from 
the ELAIS SI + Sl_5 Survey. For comparison, the black- 
shaded area represents the counts with 68% confidence con- 
tours obtained from the ISOCAM Deep/Ultra-Deep Surveys 
(Elbaz et al. 1999b). The hatched area represents the inte- 
gral counts of Serjeant et al. (2000) based on the Preliminary 
Analysis of the whole ELAIS Survey, with fluxes rescaled 
downward by an average factor of 2, as suggested by a recent 
calibration work (i.e. Vaisanen et al. 2002). The area filled 
with horizontal lines represents the range of possible expec- 
tations from no-evolution models normalized to the IRAS 
12^im local luminosity function. The asterisk represents the 
fainter end of the IRAS 12/^m source counts derived by Rush, 
Malkan and Spinoglio (1993), opportunely converted to 15 
/jm. 

Our data reduction allowed us to compute source counts 
down to fluxes ~3 times fainter than the Preliminary Analy- 
sis (PA) of the ELAIS data (Serjeant et al. 2000). In the flux 
range common to the two samples we find that, after cor- 
recting downward the PA fluxes by a factor of 2, our source 
counts are in reasonably good agreement with the counts of 
Serjeant et al. (2000) above ~2.5 mjy. However, the overall 
slope of the Serjeant et al. counts appears to be steeper than 
ours, with our counts being lower below 2.5 mjy and some- 
what higher above 10 mjy. This is probably due the fact 
that a flux dependent calibration correction, rather than a 
constant correction, should be applied to the Serjeant et al. 
(2000) fluxes. In fact, there are hints (Babbedge and Rowan- 
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Table X. The 15 /Urn Source Counts in SI + Sl_5 
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Figure 8. Differential source counts at 15 /im normalized to the Euclidean distribution (N(S) oc S~ 2,5 ). Data points (as shown also in the 
figure): A2390 (open diamonds), ISO HDF-N (open circles), ISO HDF-S (filled circles), Marano Firback (open squares), Marano Ultra- 
Deep (diagonal crosses), Marano Deep (asterisks), Lockman Deep (open triangles), Lockman Shallow (filled triangles), ELAIS SI (filled 
stars). The dotted curves are the lower and upper envelopes of our counts due to the uncertainties in the completeness curve derivation, 
as shown in figure pi The hatched area represents the range of possible expectations from no evolution models normalized to the IRAS 12 
(im Local Luminosity Function (upper limit from Rush, Malkan & Spinoglio 1993; lower limit from Fang et al. 1998). 



Robinson 2002, in preparation) that the needed correction 
factor is ~2.4 for the fainter PA fluxes and ^1.75 for the 
brighter fluxes (the factor of 2 here adopted is an average 
value) . 

Our counts, though not deep enough to detect the fast 
convergence at flux densities fainter than 0.4 mjy shown 
by the Deep/Ultra-Deep ISOCAM counts, sample very well 



the flux density range where those counts start diverging 
from no evolution models. Indeed, we observe a remark- 
able change in the slope of our counts, showing a signifi- 
cant super-Euclidean slope from about 2 mjy to ~0.45 mjy. 
Due to its large statistics, our Survey is at the moment best 
suited for determining both the exact flux density where the 
15 fim extragalactic counts steepen and the count slope itself 
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Figure 9. Integral source counts for extragalactic ISOCAM sources detected at 15 (im above a given flux S. The grey-shaded area shows the 
integral extragalactic counts with 68% confidence contours obtained from the ELAIS SI Survey (this work). For comparison, the hatched 
and the black-shaded areas represent the counts with 68% confidence contours obtained respectively from the ELAIS Preliminary Analysis, 
normalizing downwards by a factor of 2 the fluxes of Serjeant et al. (2000), and from the ISOCAM Deep/Ultra-Deep Surveys (Elbaz et 
al. 1999b). The area filled with horizontal lines represents the range of possible expectations from no-evolution models normalized to the 
IRAS 12/jm local luminosity function, as reported in figure ^. The fainter end of the IRAS 12/^m source counts derived by Rush, Malkan 
and Spinoglio (1993), converted to 15 fira, is represented by the asterisk. 



(before and after the steepening). A maximum likelihood fit 
to our extragalactic source counts with two power laws: 

dN_ ( S~ ai if S>S b , . 

dS x \ S"" 2 if S < S b l8j 

gives the following parameters: ai = 2.35±0.05, a.2 = 3.60± 
0.05, Sb =(2. 15 ±0.05) mJy. Our best fit parameters suggest 
that the steepening of the integral counts starts around 2 
mJy, then the counts keep a super- Euclidean slope down to 
the limits of our survey (^0.5 mJy). 



6 DISCUSSION 

6.1 Comparison with Deep/Ultra-Deep ISOCAM 
Surveys Source Counts 

The 15 /jm extragalactic source counts derived from the 
southern ELAIS Survey cover over two decades in flux, from 
0.5 up to 100 mJy, with a significant statistical sampling 
(325 objects). Due to the large flux density interval cov- 
ered, the ELAIS counts bridge the gap existing between the 
IRAS counts and the ISOCAM Deep/Ultra-Deep counts. 
The ELAIS Survey was planned to be a shallow survey and 



to reach an optimistic limit of about 2 mJy. As shown in Pa- 
per I, with the Lari technique we were able to go much deeper 
than expected, detecting a significant number of sources 
even at sub-mjy levels. The strength of ELAIS counts is 
at fluxes brighter than ~1 mJy, where they are highly sta- 
tistically significant and complete. At fainter fluxes, though 
the SI and Sl_5 counts are consistent over the whole flux 
range, the results are less strong due to the large incom- 
pleteness correction required by the SI data and the small 
area covered by the more complete Sl_5 data. However, the 
results are consistent with the evolution scenario found by 
other ISOCAM surveys and are able to give general hints 
on the behaviour and evolution of infrared galaxies. In par- 
ticular, ELAIS counts in the flux density range in common 
with ISOCAM Deep counts (0.5 - 4 mJy) diverge from no 
evolution models as well and steepen with a super-Euclidean 
slope (3.60 ± 0.05 for the differential form) up to the fainter 
limit. In particular, the flux density where our counts start 
diverging from no evolution predictions is (2.15 ± 0.05) mJy. 
Above this flux density, the ELAIS counts are consistent 
with no evolution, showing a slope (in differential form) of 
2.35 ± 0.05. Although similar results have been found for the 
Deep ISOCAM Surveys, between 0.5 and ~2 mJy our counts 
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are somewhat steeper and on average lower than the others. 
However, the common flux range is sampled by three Deep 
Surveys only: Marano Deep, Lockman Deep and Lockman 
Shallow. The counts drawn from the Marano Deep Survey 
appear consistent with the ELAIS counts, while the counts 
obtained from both Lockman Deep and Shallow Surveys are 
less steep at faint fluxes (especially the Lockman Shallow 
ones) and higher than ours by about a factor of 2-3 around 
2-3 mjy. 

The reason for this discrepancy is still not completely 
understood. It could be due to different separate causes or 
to a combination of them. A possible reason might be at- 
tributed to different data reduction methods applied to dif- 
ferent surveys. For example, our survey has been reduced 
with the Lari method (see Paper I), while the Marano and 
the Lockman Deep Surveys have been reduced with the 
PRETI method (Starck et al. 1999) and the Lockman Shal- 
low Survey has been reduced with the triple beam switch 
method of IAS (Desert et al. 1999). We must note that the 
PRETI method was especially designed to account for all 
the spurious effects of ISOCAM data in a more complete 
way than the triple beam switch method. A comparison be- 
tween these two methods performed in the HDF-N (Aus- 
sel et al. 1999a) has produced similar results, although not 
all the sources detected by one method were present in the 
list found by the other method. Moreover, flux densities de- 
rived for the common sources with the triple beam switch 
method were sistematically lower, by a factor of 0.82, than 
the PRETI fluxes. In fact, each method measures fluxes in 
a different way (i.e. auto-simulations for Lari, photometry 
aperture plus implicit colour correction for PRETI, fit with 
a fixed width Gaussian for the triple beam switch) and this 
might produce small differences in the photometry of the 
objects. However, if our fluxes were on the same scale as the 
PRETI ones, the assumption that the triple beam switch 
fluxes are 20% lower would furtherly increase the observed 
discrepancy between our counts and the Lockman Shallow 
counts. Viceversa, good agreement between Lockman Deep 
counts and our counts would be obtained if the Lockman 
Deep fluxes were systematically higher than ours by about 
15% (see figure 0). The flux calibration of our catalogue, as 
described in detail in Paper I, has been tested using the stars 
in the field and resulted on the same flux scale of IRAS (rel- 
ative flux scale of 1.096 ± 0.044), with photometric errors 
not larger than 10%. 

Another possible cause of the counts difference could 
be an incomplete star subtraction performed in the Lock- 
man Surveys, whose brightest flux density bins (at 2-5 mjy) 
might contain between 20 and 50% of stars (as shown in 
section ||) . 

Finally, part of the observed discrepancy could also be 
due to cosmic variance affecting small area fields, like the 
Marano Deep Field (0.2 sq. deg.) and the Lockman Deep 
and Shallow Fields (respectively 0.14 and 0.54 sq. deg.). 

6.2 Models and Interpretation 

As already mentioned, at flux densities ;J 1-2 mjy, the 
ELAIS counts are consistent with the expectation of mod- 
els assuming no evolution for extragalactic sources, while 
they strongly depart from no evolution predictions at fainter 
fluxes. The almost flat differential counts (normalised to Eu- 




T3 



S 15(im [mJy] 

Figure 10. Zoom of the normalized differential source counts at 
15 (im. The symbols are the same as in figure filled stars for 
ELAIS counts and open triangles for Lockman Deep counts. The 
upper Lockman points are the original counts, while the lower 
ones are the counts obtained with 15% lowered fluxes. 



clidean) extending from the IRAS fluxes to 1-2 mjy, followed 
by the sudden upturn below, seem to require strong evolu- 
tion for a single population rather than for the whole popu- 
lation of 15 fj,m galaxies. 

Due to the uncertainties existing in our counts at faint 
fluxes, in this paper we do not pretend to construct an evo- 
lutionary model based on our result, however, in order to in- 
terpret our data and the evolution they seem to require, we 
have compared our counts with recent evolution models for 
mid-infrared galaxies found in literature. In particular, we 
have compared our counts with the models of Xu (2000) and 
Franceschini et al. (2001). Neither model is able to reproduce 
the sharp departure feature from no-evolution predictions, 
or the low 'plateau' between 1-2 and 100 mjy shown by our 
data. The Xu (2000) model is able to fit the Deep/Ultra- 
Deep ISOCAM counts by considering a rather extreme lumi- 
nosity evolution (i.e. L(z) — L(0) x (1 + z) 4,5 ) for the whole 
infrared population. However, it largely over-predicts our 
counts above 0.8-1 mjy and its departure from no evolution 
predictions is far too smooth to reproduce the sharp upturn 
shown by our data around 2 mjy. The model predictions of 
Franceschini et al. (2001) are somewhat steeper than the Xu 
ones, but still overestimating and smoother than our counts, 
though considering a combination of luminosity and density 
evolution for star-forming galaxies only. The local luminos- 
ity functions (LLFs) on which these models are based are 
different: the one considered by Xu (2000) has been derived 
using the bivariate (15 /im vs. 60 /im luminosities) method, 
from an IRAS sample selected at 60 /im and observed by 
ISOCAM at 15 /jm, while the one used by Franceschini et 
al. (2001) is an adapted combination of the 12 fim LLF from 
Fang et al. (1998) and the bivariate 15 - 60 /im, converted to 
12 /im, from Xu et al. (1998). Franceschini et al. (2001) have 
also tried to disentangle the contributions of different popu- 
lations. Due to its larger flexibility in allowing to play with 
the different populations and their evolutionary properties, 
we have based our analysis on the Franceschini et al. (2001) 
models, trying to find a good fit to our data by varying the 
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Figure 11. Best fit model (re-adapted from Franceschini et al. 2001) to ELAIS SI source counts. The dotted line corresponds to the 
expected counts for a population of non-evolving spirals. The short-dashed line is the modeled contribution of a population of strongly 
evolving starburst (plus Seyfert 2) galaxies, while the long-dashed line is the contribution of type I AGN. The solid line represents the 
expected total source counts. 



LLF free parameters. 

These models are able to reproduce the Deep/Ultra-Deep 
ISOCAM counts by considering different evolutionary prop- 
erties for three different populations: non-evolving normal 
spirals, strongly evolving starburst plus Seyfert 2 galaxies 
and evolving AGN 1. The latter are assumed to evolve 
in luminosity as L(z) = L(0) x (1 + z) 3 up to z = 1.5 
and constant luminosity density at higher redshift. For 
the population of star-forming and Seyfert 2 galaxies, in 
a H = 50 km s _1 Mpc"\ tt m = 0.3, fi A = 0.7 Uni- 
verse, the best-fit to Deep/Ultra-Deep ISOCAM counts is 
found by Franceschini et al. (2001) by considering luminos- 
ity evolution L(z) = L(0) x (1 + z) 3,8 and density evolution 
p(L[z],z) = po(L) x (1 + z) 4 up to z break = 0.8, and no 
additional evolution at z > Zbreak- 

We have not changed either the evolution (form and rate) for 
AGN 1, nor the evolutionary scheme for star- forming (plus 
Seyfert 2) galaxies. However, for the latter, we have varied 
the luminosity and density evolution rates (respectively cv_l 
and od hereafter), zt, rea k, and the LLF normalization. The 
crucial point is the star-forming galaxy LLF and in particu- 
lar its shape at bright luminosities. In order to fit the sharp 
rise of our counts below 2 mjy, we have found essential to in- 
troduce a luminosity break (Lbre.uk) in the star-forming LLF, 



above which it quickly drops to zero. This is also strongly re- 
quired by the redshift distribution for bright sources (above 
few mjy). In fact, without any luminosity break in the LLF, 
the redshift distribution predicted by the Franceschini et al. 
(2001) model for example for sources brighter than 5 mjy, 
shows a significant secondary peak around z ~ 1, in addi- 
tion to a low redshift peak. This is not consistent with the 
redshifts measured for bright ELAIS sources (S > 2 mjy; 
Gruppioni et al. 2001; La Franca, Gruppioni, Matute et al. 
2002, in preparation), which are all found to be at rather 
low redshifts (z < 0.3 - 0.4). 

We have taken into account the redshift distribution 
constraints when looking for the best fit to our observed 
source counts. In particular, we have asked the model results 
to roughly agree with the following observational evidences: 

1. absence of high redshift peak for bright (> 2 mjy) 
sources (Gruppioni et al. 2001; La Franca, Gruppioni, 
Matute et al. 2002, in preparation); 

2. majority of sources at moderate redshifts (z < 0.3 — 0.4) 
even at fluxes S > 0.8 — 1 mjy, with a fraction of high z 
sources not larger than 30 - 35% (Pozzi, Ciliegi, Gruppioni 
et al., 2002, in preparation); 

3. redshift distribution for Deep surveys (S > 0.1 — 0.2 
mjy) showing a peak between z — 0.5 and z = 1.2 (HDF- 
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North: Aussel et al. 1999b, Aussel et al. in preparation, as 
reported by Franceschini et al. 2001; CFRS 11415+52: Flo- 
res et al. 1999). 

The best solution was found for oll = 3.0, old = 3.5, 
ztreak = 1.1, a LLF of evolving starburst population 40% 
higher than the Franceschini one and L break = 1O 1O - 8 L . In 
figure [n] the fit to source counts obtained with the above pa- 
rameters is shown, while the redshift distributions expected 
for our survey at Sis /am > 0.1, 1.0 and 2.0 mjy are plotted 
in figure These distributions are in good agreement with 
the preliminary results of optical identification for the SI 
sources on the DSS2. We find, in fact, that, while above 2 
mjy most sources have an optical counterpart brighter than 
R ~ 21.5 (92 % at Sis M m > 2 mjy and 95% at Sib y.m > 3 
mjy), between 1 and 2 mjy there is a quick drop in the 
identification fraction (it goes down to ~70%), probably as- 
sociated to a change in population (i.e. appearance of the 
high-z excess number of sources in the redshift distribution). 
At S15 > 1 mjy, the expected high-z (z > 0.6) fraction of 
sources is ~35%, in good agreement with the fraction of SI 
sources to the same flux density without optical counterpart 
on the DSS2. 

Our result shows that significant evolution is needed for 
at least a class of extragalactic objects, in order to explain 
the observed source counts below a few mjy. Above about 
10 mjy the counts are dominated by a non-evolving popula- 
tion of normal spiral galaxies, while below this flux density a 
population of strongly evolving starburst galaxies shows up 
and, rapidly rising, starts dominating the counts. At fluxes 
< 3 — 5 mjy, evolving starburst galaxies make up most of 
the observed counts, being responsible for the peak around 
0.3-0.4 mjy revealed by the Deep and Ultra-Deep ISOCAM 
Surveys. The evolution required for this class of objects is 
lower than found by Franceschini et al. (2001): oll is 3.0 in- 
stead of 3.8 and old is 3.5 instead of 4.0. However, a turnover 
at higher z (zbreak = 1.1 instead of 0.8) is needed in order 
to reproduce both the sharp rise of our counts and the faint 
flux peak of the Deep ISOCAM Surveys counts. This value 
for Zbreak is intermediate between 0.8 found by Franceschini 
et al. (2001) and 1.5 found by Xu (2000), though the lat- 
ter obtains a reasonably good fit to the Deep/Ultra-Deep 
counts and to the redshift distribution of ISOCAM sources 
in the HDF-North (Aussel et al. 1999b), by considering pure 
luminosity evolution (with oll = 4.5) for the whole 15 /im 
extragalactic population. However, the Xu (2000) model is 
not able to fit the new redshift distribution observed in the 
HDF-North with > 90% complete spectroscopic identifica- 
tion, as derived by Aussel et al. (2002 in preparation) and 
reported by Franceschini et al. (2001). 

We have obtained an estimate of the 15 /im CIRB flux 
by directly integrating the best-fit model counts down to 
Sis M m = 50 fJ-3y: 2.2 nW m~ 2 sr^ 1 . This value is in good 
agreement with the computation done by Elbaz et al. (2001, 
as reported by Chary & Elbaz 2001), who finds a value of 
2.4 ± 0.5 nW m -2 sr _1 by integrating the observed ISO- 
CAM Deep/Ultra-Deep counts down to the same flux den- 
sity limit. Our estimate of the 15 /jm CIRB flux corresponds 
to about 67% of the total resolved CIRB at 15 fim derived 
by Biviano et al. (2000) as 3.3 ± 1.3 nW m" 2 sr" 1 . 
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Figure 12. Redshift distributions predicted by our best fit model 
for 15 /im sources at different flux levels and in different areas. 
Top: S > 2 mjy, area = 4 deg 2 (SI); middle: S > 1 mjy, area = 
4 deg 2 (SI); bottom: S > 0.12 mjy, area = 6 X 10" 3 deg 2 (HDF- 
N). The latter is computed for the HDF-N area and to the same 
flux limit for a direct comparison with the redshift distribution 
observed in that field (Aussel et al. 2002, in preparation). 
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in agreement with the observed redshift distributions at dif- 
ferent flux levels (from 10 down to 0.1 mjy), predicting a 
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down to ~1.5 mjy and a rapid rise of a high-z {z me d — 1) 
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lation totally dominates below ~ 0.5 mjy. 
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